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For stars with 0.8 Mq < M < 8.0 Mq, nucleosynthesis enters its final phase during the asymptotic 
giant branch (AGB) stage. During this evolutionary period, grain condensation occurs in the stellar 
atmosphere, and the star experiences significant mass loss. The production of presolar grains can 
often be attributed to this unique stellar environment. A subset of presolar oxide grains features 
dramatic 18 O depletion that cannot be explained by the standard AGB star burning stages and 
dredge- up models. An extra mixing process, referred to as cool bottom processing (CBP), was 
proposed for low-mass AGB stars. The ls O depletion observed within certain stellar environments 
and within presolar grain samples may result from the 18 O + p processes during CBP. We report 
here on a study of the 18 0(p,7) 19 F reaction at low energies. Based on our new results, we found 
that the resonance at E^ b = 95 keV has a negligible affect on the reaction rate at the temperatures 
associated with CBP. We also determined that the direct capture S-factor is almost a factor of 2 
lower than the previously recommended value at low energies. An improved thermonuclear reaction 
rate for 18 0(p,7) 19 F is presented. 



I. INTRODUCTION 

An overwhelming majority of the matter within our 
solar system has a unique 18 0/ 16 isotopic signature. 
However, a collection of presolar grain samples feature 
peculiar oxygen isotopic ratios. These outliers are found 
within the trove of presolar grains gathered over the years 
from primitive meteorites and interplanetary dust parti- 
cles. This study is motivated by observations of preso- 
lar grains that nucleated in the atmospheres of distant, 
evolved stars before the formation of the Sun. These 
grains retain the isotopic ratios of the stellar surface they 
originated from. During the birth of the Sun, most preso- 
lar grains were annihilated as gas and dust collapsed to 
form the nascent star. As the solar system cooled and the 
Sun ascended the main sequence, the presolar grains that 
survived were incorporated into primitive meteorites. 
The study of their abnormal isotopic ratios provides cru- 
cial constraints for astrophysical models. This paper fo- 
cuses on oxide grains referred to as Group 2 grains, ap- 
proximately 15% of all presolar oxides [1]. They exhibit 
a characteristic 18 0/ 16 abundance ratio <1.5xl0~ 3 [2], 
reflecting a substantial 18 depletion [3] with respect to 
the solar value, (2.09±g;J|)x 1(T 3 [4]. 

It has been hypothesized that asymptotic giant branch 
(AGB) stars are an 18 depletion site [1]. During the 
AGB stage — the final phase of nucleosynthesis during the 
evolution of a 0.8—8.0 M star [5, 6] — a star undergoes 
substantial nucleosynthesis and mass loss. Peeling away 
the surface layers enveloping an AGB star reveals numer- 
ous burning sites and a complex interplay between these 
regions. A stellar core, composed of electron degenerate 
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carbon and oxygen, is surrounded by alternately burning 
helium and hydrogen shells. During periods of helium- 
burning, referred to as thermal pulses, thermonuclear 
runaway (TNR) occurs and drives convection between 
the two burning sites. When the TNR subsides, the star 
compensates for this period of activity by expanding and 
cooling. The hydrogen burning shell is quenched during 
expansion, and the convective envelope dredges the prod- 
ucts of nucleosynthesis to the surface of the star (third 
dredge- up). After this dredge- up event, the star con- 
tracts, and the hydrogen shell reignites. This interplay 
between the helium and hydrogen shells repeats episodi- 
cally [5]. 

During the AGB phase, 18 depletion may occur due 
to cool bottom processing (CBP) [1]. This extra mixing 
process was proposed by Wasserburg et al. [7] to account 
for isotopic anomalies, including 18 depletion, in preso- 
lar grains. During CBP, material circulates between the 
convective envelope and the radiative zone that separates 
the envelope from the hydrogen burning shell. The base 
of the convective envelope remains cool, thus distinguish- 
ing this process from hot bottom burning that occurs in 
4—7 Mq AGB stars [1, 7]. As the circulated matter ap- 
proaches the hydrogen shell, it reaches temperatures high 
enough to destroy 18 via hydrogen burning. The pro- 
cessed material is then recirculated into the convective 
envelope and transported to the stellar surface. Grains 
nucleate in the stellar atmosphere depleted in 18 O due to 
processes that occurred deep within the star. Then, pow- 
erful stellar winds inject these grains into the interstellar 
medium. The mechanism driving CBP is not understood, 
and several explanations have been proposed, including 
magnetic buoyancy [8], gravity waves [9], shear instabil- 
ity [10, 11], meridional circulation [12], and convective 
overshoot [13]. 

The observed ls O depletion in some presolar oxide 
grains and AGB stellar atmospheres helped motivate the 
introduction of CBP into AGB stellar models. These 
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models provided some insight into the class of AGB stars 
that might experience CBP and the temperature of the 
stellar plasma at the site of this extra mixing. According 
to Palmerini et al. [2], 18 O depletion by CBP may oc- 
cur in 1.5—1.7 M Q AGB stars. Cool bottom processing 
models resulted in stellar plasma temperature regimes 
within low-mass AGB stars that could allow for 18 O + p 
reactions to occur. According to Nollctt et al. [14], a 
temperature in excess of 39.8—44.7 MK (depending on 
the evolution of the core, hydrogen burning shell, and 
convective envelope boundaries) is sufficient for 18 O de- 
pletion during CBP in a 1.5 M© model star. 

The depletion of 18 O in a stellar plasma at low tem- 
peratures is driven by 18 0(p,a) 15 N and, to a lesser ex- 
tent, 18 0(p,7) 19 F. The former reaction was recently stud- 
ied indirectly by La Cognata et al. [15]. In the present 
work, we report a direct, low-energy measurement of the 
ls O(p,j) 19 F reaction. The goal of this measurement was 
to improve our knowledge of levels in the 19 F compound 
nucleus that are relevant to nuclear astrophysics. 
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FIG. 1. Truncated F level dia and 18 O + p reso- 

nances [16-21] through E]f = 844 keV. The E& b = 95 keV 
resonance corresponds to the E x = 8084 keV excited state. 
For an explanation of our reported spin and parity for this 
level, see Sec. I. Dashed arrows indicate unobserved reso- 
nances and solid arrows indicate observed resonances with 
known 7-ray decays. The proton threshold, Q P7 = 7994 keV, 
was taken from Ref. [22]. 



Within the CBP temperature regime, the 18 0(p,7) 19 F 
reaction rate may be influenced by an unobserved, low- 
energy resonance at Ejf> = 95 ± 3 keV [20, 22] (see 
Fig. 1). In the competing 18 0(p,a) 15 N reaction, a 
strength of W7 PQ = (1.6 ± 0.5) xlfT 7 eV at E}f = 95 keV 
was directly measured by Lorenz-Wirzba et al. [17]. The 
resonance strength is the integral of the reaction cross 
section and is defined as [23]: 
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where J is the compound nucleus spin, J p = 1/2 is 
the proton spin, Jt = is the 18 O target nucleus 
spin, T p is the proton partial width, T 7 is the 7- 
ray partial width, T a is the a-particle partial width, 
T x = T 7 for photon emission, and r x = r a for 
a-particle emission. In the 18 0(p,7) 19 F reaction, the 
E^ b = 95 keV resonance has never been observed, and 
none of the 7-ray decays from the resonance level are 
known. Upper limits were placed on the resonance 
strength in the past, first by Wiescher et al. [18] with 
W7 P7 < 5xl0 -8 eV and then by Vogelaar et al. [19] with 
W7 P7 < 4xl0~ 8 eV. With a proton separation energy of 
Q P7 = 7993.5994 ± 0.0011 keV [22], the Ej£ b = 95 keV 
resonance corresponds to the E x = 8084 ± 3 keV [20] 
level in the 19 F nucleus. 

The previous experimental information regarding the 
structure of this compound nucleus level is summarized in 
Table I. From the 18 0( 3 Hc,e?) 19 F experiment performed 
by Schmidt and Duhm [24], it is clear that the proton 
angular momentum transfer for the 8084 keV level is re- 
stricted to £ p = (2, 3). In La Cognata et al. [25], the 
Trojan Horse Method was used to investigate this level 
with the 2 H( 18 0,a 15 N)n reaction. They determined that 
J = 3/2 and l a = 1. Consequently, based on the angu- 
lar momentum coupling rules, the spin-parity and orbital 
angular momentum amount to 3 n = (3/2)+ and £ p = 2, 
respectively. Note that an incorrect spin and parity of 
.V = (3/2) _ was assumed previously for this level [26]. 

Here we report on a new search for the E^ b = 95 keV 
resonance in 18 0(p,7) 19 F with significantly improved sen- 
sitivity compared to previous studies [18, 19]. In the fol- 
lowing we will discuss the experimental setup, including 



TABLE I. E x 



keV level parameters. 



Parameter 


Value (eV) 


Reference 


LOJ pa 


(1.6 ± 0.5)xlO _Y 


[17] 


l^7p 7 ,UL 


< 4.0xl0" 8 


[19] 




> 1.3xl0 -11 


Sec. IIIC 


r 7 


(6.0 ± 2.5)xl0~ 1 a 


[18] 


r 


< 3.0xl0 3 b 


[17] 



a Private communication from K. Allen quoted in Wiescher et al. 
[18]. 

b Total width determined from slope of front edge of thick-target 
yield curve. 
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a brief outline of the accelerators (Sec. II A) and detector 
system (Sec. II C). Our oxygen target fabrication process 
is described in Sec. II B. The methodology we employed 
to characterize our detector efficiencies is discussed in 
Sec. IIC, and Sec. II D outlines some of the features of 
our data acquisition electronics. Results for resonant 
and non-resonant proton capture on 18 O are presented 
in Sec. Ill A and Sec. IIIB, respectively. An improved 
18 0(p,7) 19 F reaction rate is presented in Sec. IIIC. Con- 
cluding remarks are presented in Sec. IV. 



II. EXPERIMENT 
A. Accelerators 

The Laboratory for Experimental Nuclear Astro- 
physics (LENA) is dedicated to the measurement of low- 
energy nuclear reactions relevant to stellar nucleosynthe- 
sis. The cross sections measured at LENA lie within 
an energy regime that is susceptible to Coulomb sup- 
pression, and the LENA facility features key tools that 
increase the detection sensitivity. 

The LENA facility is a two-accelerator laboratory and 
consists of a high-current, low-energy Electron Cyclotron 
Resonance Ion Source (ECRIS) and an upgraded HVEC 
1MV JN Van de Graaff. The LENA ECRIS produces 
average beam currents of I p = 1.5 mA on target within 
a bombarding energy range of 50 keV < E p ab < 215 keV. 
The high-current allows for a substantial increase in low- 
energy nuclear reaction yields. The LENA 1MV JN Van 
de Graaff is capable of producing H + beam currents of 
I p < 250 /J.A at the target. Typical beam energy resolu- 
tion achieved with the JN ranges between 1—2 keV. In 
this study, it was primarily used to test our targets by 
measuring excitation functions (7-ray yield vs. bombard- 
ing energy) of the well-known 18 0(p,7) 19 F resonance at 
E^ ab = 150.82 ± 0.09 keV [21]. These excitation func- 
tions provided information on target thickness and sta- 
bility during the experiment. See Sec. II B for more in- 
formation on our 18 O targets. A detailed description of 
the LENA accelerator facility can be found in Cesaratto 
et al. [27]. 



B. Targets 

The anodic oxidation of tantalum targets was first out- 
lined by Amscl [28] . The anodization process allows tar- 
get thicknesses to be consistently reproduced, and it also 
allows the production of robust oxygen targets that re- 
main stable when exposed to intense H + beam. A new 
anodization chamber was designed and assembled for our 
measurement according to the description in Rcf. [29]. 

During fabrication, all tantalum backings were etched 
in an acid bath in order to reduce beam induced back- 
grounds by removing surface contaminants ( n B and 
19 F). Subsequently, all etched tantalum backings were 



rcsistivcly heated. These outgassed target backings were 
anodized at 64 V using 99.3% enriched 18 O water to pro- 
duce Ta 2 18 05 targets with an expected target thickness 
of -18 keV at Ej? b = 151 keV. 

Excitation functions were collected during this experi- 
ment at the well-known F}^ h = 151 keV resonance in the 
18 0(p,7) 19 F reaction with the JN Van de Graaff. Target 
thicknesses near 100 keV were estimated with the rela- 
tionship [23]: 



AE(151) _ AE(Ep) 



6eff(151) 



eff (E p ) 



(2) 



where AE is the measured target thickness in energy 
units, and e e g is the effective stopping power in the 
center-of-mass system, derived from Bragg's rule [23, 30]: 
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where M p and Misq are the mass of the proton 
and the 18 atom, eis and ex a are the laboratory 
stopping powers of protons in 18 and Ta (calcu- 
lated with SRIM [31]), and Nj are number densities 
(N = Ni6 + Nir + Ni8 Q ). We found that our tar- 
gets could withstand proton accumulations of Q p > 45 C 
without significant degradation at i^crtis _ Q.5— 1.0 mA. 



C. Detectors 

Almost all 18 0(p,7) 19 F resonances are known to de- 
cay via emission of multiple, coincident 7-rays. There- 
fore, the simultaneous detection of two or more photons 
allows an opportunity to increase the signal-to-noise ra- 
tio significantly. To accomplish this signal optimization, 
a 7-ray spectrometer consisting of several detectors was 
used. The LENA 77-coincidcnce detector system was 
assembled with an 135% HPGe detector at degrees to 
the beam and in close running geometry with the tar- 
get chamber. The distance between the HPGe detec- 
tor and the target midpoint was 1.1 cm [32]. The tar- 
get chamber and HPGe detector were surrounded by a 
16-segment Nal(Tl) annulus. Plastic scintillator paddles 
covered the two detectors on five sides and suppressed 
cosmic-ray muon events. In a two-dimensional Nal(Tl) 
vs. HPGe energy spectrum, appropriate gates were set 
during off-line data sorting. The low-energy thresholds 
set on these gates removed events caused by environmen- 
tal background ( 40 K, 208 T1), and the high-energy thresh- 
olds excluded events with a total energy that exceeded 
the excitation energy of the decaying 19 F compound nu- 
cleus. Most of the latter events were presumably caused 
by cosmic ray interactions. The LENA 77-coincidcncc 
detector was described in detail by Longland et al. [32]. 

The internal geometry of the HPGe detector is well 
known and was measured previously [33] using computed 
tomography (CT). Based on the known dimensions, rel- 
ative peak efficiencies were simulated using Geant4 [34, 
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35] by assuming mono-energetic 7-rays (E 7 = 0.05—15.0 
MeV) emitted from an extended beamspot (1.2 cm di- 
ameter) on the target. The sum-peak method [23, 36] 
was used to obtain absolute peak and total efficiencies 
for 60 Co (E 7 = 1173 keV, 1332 keV). The absolute peak 
efficiency of the LENA HPGe detector was determined 
to be Tjfsil = 0-040 ± 0.003, and the absolute total ef- 
ficiency was ij^J = 0.188 ± 0.012. Corrections for the 
finite beamspot size were estimated with Geant4 [32]. 

Peak efficiencies at lower energies were measured 
with 56 Co and at higher energies (up to 10 MeV) 
with nuclear reactions. These reactions included 
14 N(p, 7 ) 15 at E^ b = 278 keV [37] and 27 Al(p, 7 ) 28 Si at 
E^ b = 406 keV [38, 39]. All efficiency data were corrected 
for coincidence summing effects using the matrix method 
outlined by Semkow et al. [40]. The separate data sets 
were bootstrapped together across the full energy range 
to determine the experimental peak efficiencies of the de- 
tector. Between measured energies, the peak efficiencies 
were found by interpolation using Geant4. Total effi- 
ciencies were also simulated in Geant4 and then normal- 
ized to the measured 60 Co sum-peak total efficiency. The 
estimation of the 7-ray coincidence efficiencies needed for 
the data analysis is detailed in Sec. Ill A. 



D. Data Acquisition and Procedure 

As discussed in Sec. II C, the LENA spectrometer is 
composed of three different detectors: a HPGe detec- 
tor, a 16-segment Nal(Tl) annulus, and plastic scintilla- 
tor plates. Timing and energy signals were processed us- 
ing standard NIM and VME modules. The HPGe signals 
served as master triggers for the electronics. Coincidence 
and anti-coincidence events were sorted using the acqui- 
sition software JAM [41]; this system was convenient for 
applying the desired software energy and timing gates. 
Additional details concerning the electronics setup for the 
LENA 77-coincidence detector can be found in Longland 
et al. [32]. 

Initial excitation functions were produced with the 
JN Van dc Graaff at the well-known E^ b = 151 keV 
resonance for each target. In order to search for 
the E^ b = 95 keV resonance, a charge of 80 C was 
accumulated on-resonance at a bombarding energy of 
Ejf b = 105 keV, and 40 C were accumulated off-resonance 
85 keV. The average beam current on target 



at Ejf b 



butions. Sample spectra are displayed in Fig. 2, show- 
ing the on-resonance (ungated) singles HPGe detec- 
tor spectrum in blue and the coincidence gated spec- 
trum in red. For the latter spectrum, only those 
events in the HPGe detector are accepted that are co- 
incident with events in the Nal(Tl) counter of energy 



4.25 MeV < E 



Nal(Tl) 



< 10.0 MeV. It can be seen in 



Fig. 2 that this condition suppresses the environmental 
background by two orders of magnitude. Most 19 F lev- 
els decay by 77-cascades through the first (110 keV) ex- 
cited state, and all 19 F levels (with known decay schemes) 
de-excite through the second (197 keV) excited state. 
In Fig. 2, vertical dashed lines indicate anticipated lo- 
cations of the 7-rays originating from the de-excitation 
of the first excited state (110 keV) and second excited 
state (197 keV). Note that because of their low energy, 
the 110 keV photons would be significantly attenuated. 
No peaks were observed for these two secondary decays. 
In fact, although we achieved a considerably improved 
detection sensitivity compared to previous studies (by 
about half an order of magnitude; see below), no 7-rays 
from the decay of the E^ b = 95 keV resonance were ob- 
served in any of our singles or coincidence spectra. 

An improved upper limit on the resonance strength of 
the unobserved E^ b = 95 keV resonance was determined 
relative to the strength of the well-known resonance at 
Ejf b = 151 keV. The resonance strength is given by [23, 
42]: 



CJ7 = 



2e c ff A/~ m ax 

A 2 Af p BT]W 



(4) 



where e e ff is the effective stopping power at the reso- 
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III. RESULTS 

A. Resonant Data Analysis 

Gates were constructed in JAM to produce 77- 
coincidence spectra, uncover the 7-ray decay finger- 
print of the resonance, and reduce background contri- 



FIG. 2. (Color online) HPGe singles spectrum (blue) and 
77-coincidence spectrum (red). Reduction in background 
amounts to a factor of 100. The prominent background 
peak at 511 keV arises from the annihilation of pair-produced 
positrons. Dashed lines indicate the anticipated locations of 
the 1 ^ (110 keV) and 2 -> (197 keV) transitions in 19 F. 
The spectra shown represent on-resonance data, with a total 
charge accumulation of 80 C at Ep ab = 105 keV. 
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nance energy as defined in Eq. (3), A is the de Broglie 
wavelength, where [23]: 



V 
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Mr 



M t \ 4.125 x 10" 18 



M p Mt 
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(cm 2 ), (5) 



■Mnax is the total number of detected 7-rays if the target 
is considered infinitely thick, Af p is the number of incident 
protons: 



Q 

e 



(6) 



where Q is the accumulated charge on target and e is the 
unit charge in Coulomb, B is the branching ratio, rj is the 
efficiency of the detector, and W is the angular correla- 
tion. For the ratio of resonance strengths we obtain [23]: 



^795 

W7151 



05 



(7) 



151 



In this equation, W7151 = (9.7 ± 0.5) x 10 4 eV [23] from 
the weighted mean of the resonance strengths reported 
in Refs. [18, 19, 43]. All 19 F levels decay through the 
second excited state (2 — > 0), and we chose not to ex- 
clude the possibility that the 8084 keV level decays with 
a substantial primary ground state branch. Therefore, we 
used the following expression to estimate an upper limit 
for the number of 19 F compound nuclei produced [23, 44]: 



(Mnax \ 



A/rq 

Gc,P 
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Ge,P r 
^20 J 7 



(8) 



where is the simulated intensity of the 197 keV peak 
in the coincidence spectrum, Ar is the total number of 

Gc P 

simulated reactions, and r) 20 ' is the 2 — > singles peak 
efficiency (Sec. II C). This procedure was tested at the 
E^ b = 151 keV resonance, where the simulated inten- 
sities agreed with the experimental values within uncer- 
tainty (4% for the 2— >0 decay in a rectangular coincidence 
spectrum) . 

Because the E^ b = 95 keV resonance has no known 
decay scheme or branching ratios, we first calculated val- 
ues of / 7 , according to the procedure described above, 
for all bound and unbound 19 F levels with known decay 
schemes [20] . We then adopted a reasonable average from 
the ensemble of values. The statistical analysis was re- 
stricted to levels with J < 9/2 (open blue circles in Fig. 3) 
and E x > 5500 keV. The constraint on the spin was cho- 
sen to associate the calculated mean / 7 value with low- 
spin states, while the energy threshold was set so that 
the fry values were associated with complex 7-ray decay 
routes to the 19 F second excited state. As an additional 
constraint, no level with a ground state branching ra- 
tio that exceeded the total probability of emission to the 
19 F second excited state, was included in this analysis. 
This final constraint was added because the ground state 
decay mode is already included in the strength upper 
limit calculation — see Eq. (8). Results of this analysis 



for a 4.25 McV < E 



Nal(Tl) 



< 10.0 MeV gate are shown in 



Fig. 3. An average value of / 7 = 0.17 ± 0.09 represents a 
reasonable 77-coincidence correction factor estimate for 



where Aro is the upper limit on the intensity of the 
ground state transition in the singles HPGe spectrum, 
A20 is the upper limit on the intensity of the decay 
from the 19 F second excited state to the ground state 
(2 — > 0; see Fig. 1) in the coincidence-gated HPGe spec- 

Gc P 

trum, 77 RO ' is the HPGe peak efficiency for the ground 

Ge P 

state transition, r] 2Q ' is the HPGe peak efficiency of the 
2 — > transition, and / 7 is a 77-coincidence correction 
factor that depends on the 7-ray decay scheme and the 
coincidence gate selected. 

To calculate the correction factor / 7 , a Geant4 sim- 
ulation was conducted that, for a given energy level, 
used the known emission probabilities to predict the to- 
tal number of detected 7-rays arising from the 2 — > 
transition for a variety of coincidence gates. Our new 
upper limit was extracted by requiring a rectangular en- 



ergy gate of 4.25 McV < E 



Nal(Tl) 

7 



; 10.0 MeV in the 



two-dimensional Nal(Tl) vs. HPGe coincidence energy 
spectrum. The simulated coincidence histograms could 
then be sorted with the same energy gates and condi- 
tions that were used to analyze the experimental data. 
The correction factor, / 7 , was calculated by solving the 
following equation [23]: 
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FIG. 3. (Color online) 77-coincidence correction factors 
(/ 7 ) for all 19 F energy levels with J < 9/2 (open blue cir- 
cles) and E x > 5500 keV; levels with high ground state de- 
cay modes were excluded. Additionally, correction factors 
for levels with = (3/2) + are indicated by solid red cir- 
cles. These correction factors were generated by applying a 
4.25 MeV < E 7 aI(T1) < 10.0 MeV gate. The mean / 7 value for 
the entire distribution, / 7 = 0.17 ± 0.09, is represented in this 
figure by the solid blue line. The two dashed blue lines repre- 
sent the uncertainty (in this instance, the root-mean-square). 
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the E^ b = 95 keV resonance. The quoted uncertainty 
of the mean correction factor is the root-mean-square of 
the distribution. In Fig. 3, the ,V = (3/2)+ levels (lev- 
els with the same spin and parity as the Ep b = 95 keV 
resonance level) are indicated with closed red circles. 

The peak intensity upper limit for the 2— >0 transition 
(197 keV) was obtained from the HPGe coincidence spec- 
trum using the Bayesian statistical approach outlined 
in Zhu [45]. According to this method, conditional, non- 
informative posterior probability density functions were 
generated for each energy region, and peak intensity up- 
per limits were calculated. The E^ b = 95 keV resonance 
strength upper limit was then determined by generating 
normally distributed probability density functions for all 
of the other quantities that entered into the resonance 
strength calculation — Eqs. (5—8). All probability den- 
sity functions were then randomly sampled iteratively, 
and this process produced a resonance strength proba- 
bility density function that was then integrated to the 
90% confidence level. A new resonance strength upper 
limit of W795 < 7.8 xlO" 9 eV (90% CL) was obtained 
for the E^ b = 95 keV resonance in the 18 0(p,7) 19 F re- 
action. This new upper limit improves upon the upper 
limit presented in Vogelaar et al. [19] by about a factor 
of 5. 



B. Non- resonant Data Analysis 

Data collected during this study at Ep ab = 105 keV are 
also important for obtaining improved estimates for the 
direct capture cross section of 18 0(p,7) 19 F. During direct 
capture, a final bound state is created when the nucleus 
acquires a proton via 7-ray emission; this process occurs 
without the formation of a compound nucleus [46]. The 
experimental 18 0(p,7) 19 F direct capture cross section at 
glab _ 1850 keV was measured previously by Wiescher 
et al. [18]. 

While the f 1 estimation explained in Sec. Ill A re- 
lied upon a statistical argument, no such assumption 
was necessary to determine the direct capture correc- 
tion factor, /° c . To calculate the required direct cap- 
ture branching ratios, we first extrapolated the exper- 
imental cross section to Ep ab = 105 keV for all direct 
capture transitions observed in Ref. [18]. To this end, 
two different direct capture codes were employed. The 
code TEDCA [47] was used to compute the direct cap- 
ture cross section for a zero scattering potential. The 
bound state and scattering state potential parameters 
used were adopted from Iliadis and Wiescher [46]. The 
code DIRCAP [46] was utilized to perform the same calcu- 
lation with a hard-sphere scattering potential. The calcu- 
lated cross sections (from Ep m ' = 0.03—1.99 MeV) were 
normalized to the measured direct capture cross sections 
at Ejf b = 1850 keV [18]. The direct capture branch- 
ing ratios derived from this procedure were required for 
the calculation of coincidence efficiency correction fac- 
tors, /P c , using the code Geant4 (Sec. Ill A). 



No direct capture transitions were observed in any 
of the singles or coincidence spectra accumulated at 
Ejf b = 105 keV. An experimental upper limit on the total 
direct capture cross section was obtained from [23]: 



y = 



2(1 



1 



e ff Je 



E 1 ' : 



^(E) dE (10) 



-AE 



where y is the measured yield upper limit, A/20 is the 
intensity upper limit of the 2 — > transition from the 
Bayesian treatment discussed in Sec. Ill A, and ct dc (E) is 
the total direct capture cross section. This expression as- 
sumes that the effective stopping power is approximately 
constant over the target thickness, as was the case in the 
present experiment. The cross section can be rewritten 
in the form [23]: 



CT (E) = ^l e ~^v 
v ; E 



(11) 



where 5(E) is the astrophysical S-factor, E is the center- 
of-mass energy, and e -27 ™' is the Gamow factor. By as- 
suming a nearly constant S-factor over the target thick- 
ness, Eqs. (10—11) can be integrated numerically to 
extract <r(E) or 5(E) from the measured yield. This 
set of calculations was performed for the same 77- 
coincidence gate used in Sec. Ill A. For the total exper- 
imental astrophysical S-factor, we found an upper limit 
of S£gj ^ 8 -l kcV b ( 90% CL )> corresponding to a di- 
rect capture cross section upper limit of c^tai < 1.8 pb 
(90% CL). Note that these values are nearly independent 
(within 2%) of the direct capture code used to calculate 
the branching ratios at Ep ab = 105 keV. 



J2 

> 



CO 



18 0(p,Y) 19 F 

j°> Wiescher et al. (1 980) 




200 400 



600 800 1000 1200 1400 1600 1800 2000 
r-c.m. . 



(keV) 



FIG. 4. (Color online) Total direct capture S-factor for 
18 0(p,7) 19 F. The solid lines represent direct capture model 
calculations: (black) Wiescher et al. [18]; (blue) using code 
DIRCAP; (red) using code TEDCA; the latter two results are 
normalized to the measured direct capture cross section at 
Ejf b = 1850 keV [18]. Our measured upper limits (90%, 95%, 
99% confidence levels) at Ep m ' = 99.4 keV are displayed as 
three black arrows. 
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Our experimental total S-factor upper limit (90% CL) 
at Ep ab = 105 keV is shown in Fig. 4, along with the 
values corresponding to the 95% and 99% confidence lev- 
els. It is interesting to compare our measured upper limit 
values with direct capture model calculations. The black 
solid curve represents the total S-factor reported by Wi- 
escher et al. [18], while the red and blue solid lines were 
calculated in the present work using the codes TEDCA [47] 
and DIRCAP [46], respectively. The latter two were nor- 
malized to the previously measured direct capture cross 
section at E]f b = 1850 keV [18]. At Ejf b = 105 keV, 
our measured upper limits are smaller than the predic- 
tion of Wiescher et al. [18] by about a factor of 2. The 
DIRCAP S-factor (blue line) was only marginally consis- 
tent with our experimental upper limit (90% CL) while 
the extrapolation derived from the code TEDCA (red line) 
fell within the 90% confidence level. 



C. Reaction Rates 

Thermonuclear reaction rates for 18 0(p,7) 19 F were 
calculated using the Monte Carlo method of Longland 
et al. [48]. In the Monte Carlo calculation, using the 
code RatesMC, we adopted the same nuclear physics in- 
put as in Ref. [49], except for the E^ b = 95 keV res- 
onance strength, the total direct capture S-factor, the 
Q- value [22], and the resonance energies. 

Based on Iliadis and Wiescher [46] , we adopted the 
TEDCA extrapolation of the total direct capture S-factor, 
normalized at Ep ab = 1850 keV [18]. For bombarding 
energies below E£' m ' = 2.0 MeV, our adopted total S- 
factor can be expanded around E = 0, with the result: 



5(E) « 5(0) + S'(0)E 
= 7.06 + 2.98 x 10 _3 E 



^"(0)E 2 



(12) 



2.60 x 10 -7 E 2 (keV b), 



where E is the center-of-mass energy. Note that at low 
energies, our new direct capture S-factor is significantly 
smaller than the result reported in Wiescher et al. [18]. 
Also note that the S'(0) coefficient presented in Wiescher 
et al. [18] was reported incorrectly and should in fact be 
S'(0) = -0.34 xlO -3 b [50]. This correction is already 
applied to the black line displayed in Fig. 4. 

In the present work, we reported on an improved 
upper limit of the Ep b = 95 keV resonance strength, 
W7 < 7.8 xl0~ 9 eV (90% CL). For this particular 19 F 
level, we may also estimate a lower limit for the resonance 
strength based on the available resonance properties (see 
Tab. I). The ratio of resonance strengths in the (p,7) and 
(j),a) channels, according to Eq. (1), is given by: 



(13) 



The (p,a) strength was measured by Lorenz-Wirzba et al. 
[17], with the result w~/ pa = (1.6 ± 0.5) xlO -7 cV. An 
upper limit for the total width of T < 3xl0 3 eV was 



obtained from the slope of the low-energy edge of the 
thick-target yield curve [17], implying an upper limit of 
L Q < 3xl0 3 eV for the a-particle partial width. Finally, 
a value of T 7 = (6.0 ± 2.5) xl0 _1 eV was reported for 
the 7-ray partial width in Wiescher et al. [18]. With 
these input values and their associated uncertainties, we 
found, from Eq. (13), a lower limit on the (p,7) strength 
of w 7p7 > 1.3xl0- n cV. 

Since we were able to calculate both an upper and a 
lower limit on the strength, we estimated a recommended 
value and a factor uncertainty using the following equa- 
tions [48]: 



ujj = V^7ll x cl>7 UL = 3.2 x 10 



-10 



eV, 



f.u. = 4 /^=25. 

^7LL 



(14) 
(15) 



In our Monte Carlo procedure, the rate contribution of 
the E^ b = 95 keV resonance strength was found by 
randomly sampling a lognormal distribution constructed 
from the mean value and factor uncertainty. Associat- 
ing resonance strengths with lognormal distributions is 
discussed in Longland et al. [48]. 

Our new low, median, and high 18 0(p,7) 19 F reaction 
rates (corresponding to 0.16, 0.50 and 0.84 quantiles, re- 
spectively, of the cumulative reaction rate distribution) 
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FIG. 5. (Color online) (Left) Reaction rate probability den- 
sity functions (red) for 18 0(p,7) 19 F at 0.02 GK, 0.2 GK, and 
2.0 GK generated by the RatesMC Monte Carlo code. The 
lognormal approximations are overlaid in black. (Right) The 
corresponding cumulative probability functions used to define 
the low, median and high rates as 0.16, 0.5, and 0.84 quan- 
tiles, respectively. 
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T (GK) Low rate 



0.010 
0.011 
0.012 
0.013 
0.014 
0.015 
0.016 
0.018 
0.020 
0.025 
0.030 
0.040 
0.050 
0.060 
0.070 
0.080 
0.090 
0.100 
0.110 
0.120 
0.130 
0.140 
0.150 
0.160 
0.180 
0.200 
0.250 
0.300 
0.350 
0.400 
0.450 
0.500 
0.600 
0.700 
0.800 
0.900 
1.000 
1.250 
1.500 
1.750 
2.000 
2.500 
3.000 
3.500 
4.000 
5.000 
6.000 
7.000 
8.000 
9.000 
10.000 



2.923x10" 
2.581x10" 
1.719x10" 
9.299x10" 
4.223x10" 
1.709x10" 
6.180x10" 
6.283x10" 
4.820x10" 
3.235x10" 
8.728x10" 
1.206x10" 
9.349x10" 
9.553x10 
3.573x10 
5.495x10 
4.553x10 
2.437x10 
9.505x10 
2.921x10 
7.484x10 
1.663x10 
3.300x10 
5.970x10 
1.581X10" 1 
3.388X10" 1 
1.274x10° 
2.932x10° 
5.153x10° 
7.695x10° 
1.037X10 1 
1.303X10 1 

1.841 xlO 1 
2.464X10 1 
3.356X10 1 
4.759X10 1 
6.916X10 1 
1.719xl0 2 
3.630X10 2 
6.403 xlO 2 
9.921 xlO 2 

1.842 x10 s 
2.798 x10 s 
3.777x10 s 
4.758x10 s 
6.600 x10 s 
(8.727xl0 3 ) 
(1.167xl0 4 ) 
(1.452xl0 4 ) 
(1.718xl0 4 ) 
(2.032xl0 4 ) 



Median rate 
4.967xl0 -24 
4.222X10" 23 
2.676xl0" 22 
1.370xl0~ 21 
5.973xl0" 21 
2.321x10 
8.151x10 
8.027x10 
6.045x10 
4.041x10 
1.093x10 
1.555x10 
1.227x10 
1.274x10" 
4.764xl0" 7 
7.304xl0" 6 
6.049xl0" 5 
3.239x10 
1.263x10 
3.880x10 
9.946x10 
2.211x10 
4.386x10 
7.942x10 
2.102x10 
4.507x10 
1.694x10 
3.903x10' 
6.853x10° 
1.020X10 1 
1.370X10 1 
1.715X10 1 
2.387X10 1 
3.121X10 1 
4.137X10 1 
5.709 xlO 1 
8.167X10 1 
2.000 xlO 2 
4.213X10 2 
7.430 xlO 2 
1.149x10 s 
2.129x10 s 
3.230x10 s 
4.369 x10 s 
5.507x10 s 
7.729 x10 s 
(1.056xl0 4 ) 
(1.411xl0 4 ) 
(1.757xl0 4 ) 
(2.078 xlO 4 ) 
(2.458xl0 4 ) 
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-3 
-3 
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-2 
-2 
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High rate 
8.785xlO -24 
7.344xl0 -23 
4.468xl0" 22 
2.181xl0 -21 
8.985xl0" 21 
3.309x10 
1.113x10 
1.041x10 
7.721x10 
5.168x10 
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-19 
18 
18 
16 
14 
12 
10 



1.402x10 
2.340x10 
2.056x10 
1.818xl0" 8 
6.503xl0" 7 
9.891xl0" 6 
8.137x10" 
4.347x10" 
1.692x10" 
5.191x10" 
1.330xl0" 2 
2.957xl0" 2 
5.867xl0" 2 
1.062X10" 1 
2.814X10" 1 
6.033X10" 1 
2.266x10° 
5.212x10° 
9.134x10° 
1.360X10 1 
1.819X10 1 
2.269X10 1 
3.123X10 1 
3.988X10 1 
5.129X10 1 
6.938X10 1 
9.819X10 1 
2.380xl0 2 
4.975 xlO 2 
8.726xl0 2 
1.342 x10 s 
2.487x10 s 
3.769 x10 s 
5.130x10 s 
6.507x10 s 
9.353 x10 s 
(1.277xl0 4 ) 
(1.707xl0 4 ) 
(2.125xl0 4 ) 
(2.514xl0 4 ) 
(2.974 xlO 4 ) 



are tabulated in Tab. II over a stellar plasma tempera- 
ture range of 0.01 — 10.00 GK. Reaction rate probability 
density functions at a few sample temperatures (0.02, 
0.2, 2.0 GK) are displayed as red histograms in Fig. 5 
(left panel), with the lognormal approximations shown 
as black solid lines. On the right, the corresponding cu- 
mulative probability functions are shown with the dashed 
lines indicating the low, median, and high rates. It can be 
seen that a lognormal function approximates the actual 
Monte Carlo distribution rather well. 

Figure 6 compares our new reaction rate with the one 
published by Iliadis et al. [49] . The new (solid lines) and 
previous (dotted lines) high and low rates are normalized 
to the previous recommended rate [49]. Note that the 
previous rates contained two small mistakes: (i) an erro- 
neous assignment of = (3/2)", and (ii) the incorrectly 
reported value of S'(0) = 0.34xl0" 3 b from Rcf. [18]. 
The dashed vertical line at 44.7 MK indicates the high- 
est temperature threshold at which, according to Nollett 
et al. [14], CBP can occur. The vertical dashed line at 
5.5 GK represents the stellar temperature beyond which 
the rates must be found with the aid of Hauser-Feshbach 
calculations. This threshold was computed based on the 
methodology outlined by Newton et al. [52]. 

The difference, in Fig. 6, between new and previous 
rates at temperatures below 50 MK can be explained by 
our lower estimates both for the contributions from di- 
rect capture and the resonance at E^ b = 95 keV. Since 
the new rates are smaller at CBP threshold tempera- 
tures compared to the previous results, it is even less 
likely that the 18 0(p,7) 19 F reaction contributes signifi- 
cantly to the depletion of 18 observed in stellar atmo- 
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TABLE II. Experimental Monte Carlo-based ls O(p,7) 19 F re- 
action rates (in units of cm mol" 1 s" 1 ). For T > 5.5 GK, 
rates were matched to Hauser-Feshbach results [51]. 



FIG. 6. Present (solid lines) and previous [49] (dotted lines) 
high and low reaction rates, normalized to the recommended 
previous rates. The vertical dashed line at 44.7 MK repre- 
sents the highest lower limit on CBP temperatures within 
a low-mass AGB star according to Ref. [14]. The vertical 
dashed line at 5.5 GK represents the temperature at which 
the experimental rates need to be extrapolated with the aid 
of Hauser-Feshbach results [51]. 
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spheres and presolar grain samples. The slight increase 
in the rate above 44.7 MK is dependent upon the cal- 
culated E|f b = 95 keV strength upper limit and cannot 
account for observed 18 depletions. The difference at 
temperatures in excess of 5 GK is solely caused by our 
treatment of the direct capture contribution: the S-factor 
expansion was artificially cut off at E^™- = 1.0 MeV in 
previous work [18, 49], while in the present work the S- 
factor is calculated up to energies of E^ 111 - = 2.0 MeV 
(Fig. 4), resulting in a much higher cutoff value and a 
significantly increased direct capture contribution. 

The fractional contributions to the total 18 0(p,7) 19 F 
reaction rates are shown in Fig. 7 for low temperatures. 
It can be seen that our lower estimate for the direct 
capture process (blue solid line) contributes significantly 
(>10%) at temperatures of 0.011-0.05 GK. Our consid- 
erably lower estimate for the Ep b = 95 keV resonance 
(red solid line) yields a negligible contribution (<3%) in 
the temperature region relevant to cool bottom process- 
ing. The dashed lines represent the reaction rate con- 
tributions of the other resonances (Ep b = 22, 151 and 
844 keV) that influence the rate at stellar plasma tem- 
peratures below 0.2 GK. 

The ratio of 18 0(p,a) 15 N and 18 0(p,7) 19 F high and 
low rates is shown in Fig. 8. The dotted lines are based 
on the results of Ref. [49] alone, while the solid lines incor- 
porate the new 18 0(p,7) 19 F rates. For the temperature 
region relevant to CBP, the established (p,a) rate [49] 
exceeds the (p,7) rate by a factor of 5100—1700 over the 
temperature range 0.03—0.05 GK. From our improved 
E|| b = 95 keV resonance strength upper limit and our 
refined direct capture S-factor, we support the conclusion 
that the (p,7) reaction does not contribute significantly 
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FIG. 8. Ratios between (p,a) low and high reaction rates 
from Iliadis et al. [49] and the present (p,j) high and low re- 
action rates, respectively (solid black lines). The correspond- 
ing ratios based solely on the previous rates [49] are shown as 
dotted lines. The vertical dashed line at 44.7 MK indicates 
the highest CBP temperature threshold according to Nollett 
et al. [14]. 



to the overall 18 destruction at temperatures suggested 
for CBP to occur in low-mass AGB stars. Future efforts 
to study 18 depletion by CBP in AGB stars should fo- 
cus on direct measurement of the 18 0(p,a) 15 N reaction 
at low energies. 
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FIG. 7. (Color online) Fractional contributions at low tem- 
peratures (Tg < 0.2) to the 18 0(p,7) 19 F reaction rate are 
shown. The blue line represents the direct capture contribu- 
tion, the red line is the EiJ b = 95 keV contribution, and the 
dashed, dashed-dotted and dashed-double-dotted lines are the 
contributions from the E^ b = 22 keV, E^ b = 151 keV and 
E^ b = 844 keV resonances, respectively. 



IV. CONCLUSION 

A study of the 18 0(p,7) 19 F reaction was performed 
at the Laboratory for Experimental Nuclear Astro- 
physics (LENA). A new resonance strength upper limit 
of oj-/ < 7.8 xlO" 9 cV (90% CL) for the E|f b = 95 keV 
resonance was measured that improves upon the previ- 
ous (p,7) upper limit published by Vogelaar et al. [19] 
by about half an order of magnitude. Our data also 
allow for a significant improvement of the total direct 
capture S-factor prediction. Our direct capture S-factor 
amounts to about half of the previously accepted value 
at low energies [18]. With this experimental information, 
new Monte Carlo-based reaction rates for 18 0(p,7) 19 F 
are derived. We find that the new reaction rates in the 
hypothesized CBP temperature regime are even smaller 
than previously assumed. Clearly, 18 depletion in low- 
mass AGB stellar atmospheres and some presolar oxide 
grains is dominated by the competing ls O(p,a) 15 N re- 
action. Future studies of 18 O depiction by cool bottom 
processing in low-mass AGB stars should focus on direct 
measurement of the (p,a) reaction at low energies. 
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